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Occurrence of the Berezinskii-Kosterlitz-Thouless (BKT) transition is investigated by superfluid 
density measurements for two-dimensional (2D) disordered NbN films with disorder level very close 
to a superconductor- insulator transition (SIT). Our data show a robust BKT transition even near 
this 2D disorder-tuned quantum critical point (QCP). This observation is in direct contrast with 
previous data on deeply underdoped quasi-2D cuprates near the SIT. As our NbN films approach 
the quantum critical point, the vortex core energy, an important energy scale in the BKT transition, 
scales with the superconducting gap, not with the superfluid density, as expected within the standard 
2D-XY model description of BKT physics. 

PACS numbers: 74.40.-n, 74.40.Kb, 74.62. En, 74.25.Ha 



I. INTRODUCTION 

Berezinskii-Kosterlitz-Thouless (BKT) transition, the 
only phase transition which can occur in the two- 
dimensional XY (2D-XY) model without breaking the 
continuous symmetry of the model, 1 has generated great 
interest in condensed matter community for many years. 
It has been used to describe the superconductor-to- 
normal-metal thermal phase transition in 2D supercon- 
ducting films in the context of free vortices emerging 
from a bath of thermally excited vortex-antivortex (V- 
aV) pairs, instead of breaking of Cooper pairs themselves. 

There are several predicted experimental signatures of 
this transition. 2 For example, above the transition tem- 
perature, the coherence length would diverge exponen- 
tially in the distance from the transition instead of the 
usual power-law, leading to a peculiar temperature de- 
pendence of the resistivity above the transition. 3 But this 
temperature range is usually exceedingly small. 3 ~ 5 The 
most direct and convincing evidence is that, universal 
and discontinuous drop in superfluid density is expected 
at the transition temperature. 6 This has been shown 
beautifully in superfluid helium-4 system. 7 In ultrathin 
conventional (Mo 77 Ge23, 8 InC^ 9 ' 10 , NbN 5 ' 11 ) and quasi- 
2D cuprate (YBa 2 Cu 3 07- :!; 12 ~ 14 , Bi 2 Sr 2 CaCu 2 8+:r 15 ) 
superconducting films, however, results arc rather com- 
plex: (1) While significant drops in superfluid densities 
are indeed observed, they do not occur where 2D-XY 
model predicts - in ultrathin conventional films, they 
occur earlier than expected. 5,8,9,11 (2) In strongly un- 
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derdoped cuprates, thick films of YBa 2 Cu307 



Bi 2 Sr 2 CaCu 2 08+;r, and crystals of YBa 2 Cu 3 07 
thermal critical fluctuations are not observed near T c , 
even though samples near optimal doping do exhibit crit- 



ical fluctuations. 12 ' 15 It is worth noting that within the 
context of layered cuprates the possibility to identify 
BKT features associated to each bilayer unit relies on 
the general expectation that layers are weakly coupled. 
Thus, it is particularly surprising that when cuprate films 
arc underdoped to near a superconductor-insulator tran- 
sition (SIT), any thermal critical behavior, evidenced 
by the sharp downturn of superfluid density, disappears 
and the T-dependence of superfluid density goes quasi- 
linearly with the temperature all the way to T c 13 ~ 15 . This 
contradicts the fact that underdoping usually increases 
anisotropy in cuprates, so that thermal critical behavior 
should be more robust than its counterpart near optimal 
doping. 

All the observations above indicate that there are sev- 
eral physical mechanisms at play in 2D superconducting 
films that are not captured by the 2D-XY model descrip- 
tion of the BKT physics. First, as pointed out by one 
of the authors here, the relative energy scales involved 
in the BKT transition might not be universal after all 16 . 
For instance, the experimental data can be described well 
by allowing the ratio between the vortex core energy fj, 
and the superfluid density to deviate from the 2D-XY 
model value. A smaller (compared to the 2D-XY model 
prediction) or a larger fj, can be used to fit the data of 
ultrathin conventional superconducting films 5 ' 11 and lay- 
ered cuprates 17,18 respectively, to account for the early 
or late drop in superfluid density. 

Second, as evidenced mainly by scanning tunneling mi- 
croscopy, intrinsic inhomogeneities emerge in these films 
of both conventional 19-22 and cuprate 23 superconduc- 
tors especially when they are underdoped or driven to a 
very high disorder level. These inhomogeneities tend to 
broaden the transition and smear out the discontinuous 



drop. Any quantitative analysis then must take the local 
distribution of the superfiuid densities into account 18 ' 24 . 
This will also complicate the analysis. 

Third, when the system is pushed near the verge of 
a SIT, no matter by disorder or undcrdoping, there will 
be quantum fluctuations near such a quantum critical 
point (QCP). For example, there can be quantum V-aV 
pairs existing even at zero temperature. 25 How quantum 
fluctuations affect BKT transition is still an open ques- 
tion. Goldman et al. suggested recently that macro- 
scopic quantum tunneling in non-uniform thin films pre- 
vent resistivity dropping to zero below BKT transition. 26 
Previous studies on deeply underdoped cuprates also 
show that thermal critical behavior (BKT physics) per- 
sists in ultrathin films 27 but it disappears in thick 
films near the QCP ,15 . This phenomenon has been 
shown to be universal in cuprates because it is robust 
against huge differences in anisotropy (YBa2Cu307_ x vs 
Bi2Sr2CaCu208+;c) and disorder (thick films vs crystals). 
It seems BKT physics surrenders to quantum effects near 
a quantum critical point. 

The main purpose of this article is to study the evo- 
lution of BKT physics near a QCP in 2D conventional 
superconducting films. While there are some previ- 
ous superfiuid density studies on 2D superconducting 
films, none of these films have been pushed to extremely 
high disorder so that a superconductor-to insulator (SIT) 
transition can be seen. In our case, by reducing the thick- 
ness and adding disorder in NbN films, we are able to 
drive NbN films with T c ~ 15K smoothly all the way to 
insulating. We are able to push the disorder smoothly to 
SIT and still have reasonably sharp transitions. Temper- 
ature dependences of superfiuid densities in these films 
are measured by a two-coil apparatus. Qualitatively, 
BKT transitions are observed for all the films even on 
the verge of SIT. By analyzing the data within the same 
theoretical scheme proposed in a previous article, 5 we 
also observed an increase of both vortex-core energy and 
inhomogeneities as function of disorder, that consistently 
extend the previous studies on intermediately disordered 
films to highly disordered films near a 2D-QCP. This in- 
dicates that BKT physics remains robust against high 
disorder level or other quantum effects near a QCP in 
conventional 2D superconducting films. This robustness 
is in direct contrast with similar studies on layered un- 
derdoped cuprates, where BKT physics vanishes. We 
conclude this difference is because in deeply underdoped 
layered cuprates the increase of both the vortex-core en- 
ergy and of the coherence length can conspire to mask 
the occurrence of 2D behavior near the SIT. 



II. EXPERIMENTAL 

The superconducting NbN films were deposited by re- 
active magnetron sputtering of a pure Nb target in an 
Ar+N 2 gas mixture, at a total pressure of about I0~ 3 
mbar. The epipolished R-plane sapphire substrates were 
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FIG. 1: Sheet resistance R, normalized by quantum resistance 
h/Ae 2 , vs. temperature for many NbN films with different 
thickness d and disorder. Superconductor-insulator transition 
is observed for films with R(20K) < h/Ae 2 (black curves). 
Some films (red curves) are very close to SIT but the mea- 
surement temperatures go down only to 4.2K. 



kept at 550° C during the film growth. Deposition rate 
(0.17nm/sec) is calibrated so thickness of the film is in- 
ferred from the sputtering time. The deposition process 
was optimized with respect to the partial pressure of N2 
and the deposition rate to provide the highest transition 
temperature for films with the smallest studied thick- 
ness. More details of growth can be found in this paper 28 . 
Many physical parameters have been measured for these 
NbN films. 29 They are patterned to nanowires and used 
to make single photon detectors. 28 We emphasize that 
the films are homogeneously disordered because (I) Sheet 
resistance is almost a constant above the transition tem- 
perature and does not show any discontinuity at higher 
temperature. (2) Conducting films with reasonably sharp 
and single transition can be grown with the thickness of 
only three or four unit cell. (3) It is generally easier to 
be homogeneously disordered for a binary compound, like 
well-studied InO^ and TiN films. 

Figure I shows that a nice superconductor-insulator 
transition is observed when the film thickness is reduced 
to a few unit cell 30 (one unit cell = 0.44nm). Supercon- 
ducting and insulating films are separated by quantum 
resistance for cooper pairs, which is h/Ae 2 — 6.45 kO. 
Thick NbN films have T c sw 15K, close to the bulk value. 
As films get thinner or more disordered, sheet resistance 
increases and T c drops. Near SIT, we are able to con- 
sistently reproduce films with T c about 4K with sheet 
resistance ~ 5.5/cO above T c . Unfortunately, these ultra- 
thin films degrade in the air and the effective sheet resis- 
tance will change. Both the superfiuid density and T c will 
change too. This degradation prevents us from directly 
comparing resistance and superfiuid data on the same 
film. But it gives us another way to tune the SIT with- 
out sample-to-sample variation. As we will show later, in 
terms of SIT tuning, there is no difference among thick- 



ness, disorder or degradation. 

Superfluid densities are measured by a two-coil mutual 
inductance apparatus. 31 The film is sandwiched between 
two coils, and the mutual inductance between these two 
coils is measured at a frequency cu/2ir = 50 kHz. The 
measurement actually determines the sheet conductivity, 
Y = (<7i + i<J2)d, with d being the superconducting him 
thickness and a being the conductivity. Given a mea- 
sured him thickness, a is calculated as: a = Y/d. The 
imaginary part, a%, yields the superfluid density through: 
wcr 2 = n s e 2 /m, which is proportional to the inverse pen- 
etration depth squared: A~ 2 (T) = nousa2 (T) , where po 
is the permeability of vacuum. As is customary, we refer 
to A~ 2 as the superfluid density. The dissipative part of 
the conductivity, <Ti(T), has a peak near T c , whose width 
provides an upper limit on the spatial inhomogeneity of 
T c over the 10 mm 2 area probed by the measurement. 
Data are taken continuously as the sample slowly warms 
up so as to yield the hard-to-measure absolute value of 
A~ 2 and its T-dependence. This two-coil technique is 
powerful 16 to study thermal critical behavior like BKT 
transition near T c . It is also unique for 2-D films because 
it can give sheet superfluid density d/X 2 (T) without the 
knowledge of the film thickness, which is the case here. 

Mutual inductance data, MI(T)/MI(15 K), of some 
films are shown in Figure 2. Nice sharp transitions are 
observed for most films. As films get thinner, or get more 
disordered, or simply degraded, T c drops and the normal- 
ized mutual inductance at T -C T c grows. This shows 
that the ability of film to screen magnetic field gets re- 
duced as the film get thinner. This ability is directly 
related to the areal superfluid density of the film. The 
imaginary part of the mutual inductance, which shows a 
dip at the transition, relates to the dissipation of moving 
vortices. The width of the dip shows that inhomogeneity 
gets larger as the film disorder is increased to near SIT, 
consistent with other local gap measurements. 20 



III. THEORETICAL ANALYSIS OF THE DATA 

We start by briefly reviewing the basic elements of the 
BKT transition that will be needed for the fitting of the 
experimental data. 16 As we mentioned at the beginning, 
the BKT transition was originally formulated within the 
context of the two-dimensional (2D) Ay-model, which 
describes the exchange interaction between classical two- 
component spins with hxed length 5=1: 
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where J is the spin-spin coupling constant and 9i is the 
angle that the i-th spin forms with a given direction, and 
i are the sites of a square lattice. Within the context of 
2D superconductors 8i plays the role of the SC phase, and 
J (now written as J s , and referred to as the superfluid 
stiffness) is connected to the areal density of superfluid 
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FIG. 2: Temperature dependence of mutual inductance, 
MI(T)/MI(15 K), including both real (upper curves) and 
imaginary parts (lower dips), of many NbN films with dif- 
ferent thickness and disorder level. For clarity purposes, not 
all data are shown. 
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Usually both quasiparticle excitations and phase fluctu- 
ations contribute to the depletion of J s towards zero. In 
the case of our NbN films the quasiparticle contribution 
can be well accounted by the dirty-limit BCS expression, 



J BCb {T) _ A(T) 
J BCS {0) ~ A(0T 



tanh 



A(T) 
2k B T 



(3) 



by using eventually A(0)/Tbcs as a free parameter, to 
account for the relatively large A(0)/T c ratio reported in 
NbN as disorder increases. 20 . 

For what concerns transverse (i.e. vortical) phase fluc- 
tuations their effect will be accounted for by numer- 
ical solution of the BKT renormalization-group (RG) 
equations, whose relevant variables are the dimensionless 
quantities 1 ' 2 ' 16 : 



n,J BCS (T) 
T 
g(0) = 2ne-^, 



K(0) 



(4) 
(5) 



where /i is the free energy of a vortex core, with radius 
about equal to the superconducting coherence length, 
£(T), and g is called the vortex fugacity. Notice that 
jBCS en t ers here to determine the initial value of K, i.e. 
its short-distance value. Its long distance value follows 
by the solution of the well-known RG equations 1 ' 2 ' 16 : 
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FIG. 3: Role of the vortex-core energy on the BKT transition. 
The solid black line represents the temperature dependence of 
1/A 2 within BCS theory, as described by Eq. (3), for typical 
parameter values appropriate for NbN films. The BKT tran- 
sition temperature depends on the value of the vortex-core 
energy. For fi as large as in the XY model the transition oc- 
curs (dotted green line) at the intersection between the BCS 
curve and the universal line 2T/n. However, for smaller [i val- 
ues, J s is renormalized with respect to its BCS counterpart 
already before the transition, so that the transition occurs at 
a lower Tbkt (see dashed blue line and red dot-dashed line). 
Notice that in all these cases the universal relation (9) is satis- 
fied, and 1/A 2 jumps discontinuously to zero after intersection 
with the 2T/% line. 



where I = ln(r/£) is the rescaled length scale. The 
observed superfluid density is identified by the limiting 
value of K as one goes to large distances 6 : 
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The basic idea of the RG equations is to look at the 
large-scale behavior of the superfluid stiffness and of the 
vortex fugacity. When g — > it means that single-vortex 
excitations are ruled out from the system, which is then 
SC: indeed, as one can see from Eqs. (6)-(7) when g — » 0, 
K goes to a constant and then J s from Eq. (8) is finite. If 
instead g — > oo at large distances it means that vortices 
proliferate and drive the transition to the non-SC state, 
since K — > 0. The large-scale behavior depends on the 
initial values of the coupling constants K, g, which in 
turn depend on the temperature. The BKT transition 
temperature is defined as the highest value of T such 
that K flows to a finite value, so that J s is finite. This 
occurs at the fixed point K = 2, g = 0, so that at the 
transition one always has: 



va l ~, t \ o . kJ 8 {T B kt) ., 
K (I -> oo, ±bkt) = 2, =>■ — — = 2, 
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while above it, J s = 0. As a consequence, at TbkTi Js 
jumps discontinuously from the universal value ITbkt /tt 



FIG. 4: Role of inhomogeneity on the BKT transition. The 
experimental data (black curve) correspond to the sample la- 
beled as S193vl in Table I. While the median J (T) (blue 
dot-dashed line) of the Gaussian distribution (11) of possi- 
ble Ji realizations has a sharp transition, the average stiffness 
Jav{T) from Eq. (10) vanishes with a smoother tail. The de- 
viation from the BCS curve (red dashed line) before than the 
transition is due instead to the low value of the vortex-core 
energy, see Fig. (3). 



to zero. However, it should be emphasized that al- 
ready before Tbkt the effect of short length-scale vortex- 
antivortex pairs is in general to deplete J s with respect 
to its initial value, given by the BCS estimate (3). This 
effect is usually negligible when /i is large, as it is the case 
for superfluid films 7 or within the standard XY model 16 , 
where jAxy ~ (tt 2 /2) J s . In this case one can safely esti- 
mate Tbkt as the temperature where the line 2T/ir in- 
tersects the J BCS (T) from Eq. (3), see Fig. 3. However, 
as [i decreases the renormalization of J s due to bound 
vortex pairs increases, and consequently the deviation of 
J s from its BCS counterpart starts considerably before 
the transition temperature itself, 16,17 see Fig. 3. 

At intermediate and strong disorder STS experiments 
have shown that NbN films 20-22 exhibit a spatial inhomo- 
geneity of the SC spectra, that becomes particularly pro- 
nounced near the SIT. Even though STS spectra probe 
only the local DOS of the sample, one would expect that 
the same inhomogeneity reflects also in the local super- 
fluid stiffness. As a consequence, one can imagine that 
the sample admits a given distribution of local Jj values 
with probability density P{Ji) and local BCS and BKT 
transition temperatures T l c and T 1 B kt^ respectively. A 
possible phenomenological way to estimate the overall 
superfluid stiffness 24 is to compute the average J av as: 



j av (T) = j2p(wi(n 



(10) 



where P{Ji) can be taken for example as a Gaussian dis- 
tribution centered around the experimental value of Jq 



at T = 0, 



P(Ji) 



27TCT 



cxp[-(Ji- Jo) 2 /2a 2 ] 



(11) 



When Ji = J the corresponding J l 8 {T) = Jq(T) coin- 
cides with the BKT curve obtained from the BCS fit 
(3) of the experimental data, shown with a dot-dashed 
line in Fig. 4. For the remaining Ji values we rescale 
the corresponding BCS temperatures as Ji/T BCS ~ 
Jbcs(0)/Tbcs and we compute J l s (T) and the corre- 
sponding BKT temperature T BKT by the numerical so- 
lution of the RG equations (4)-(5) above. Once obtained 
this set of JI (T) curves we compute at each temperature 
the average value J av (T) according to Eq. (10). When all 
the stiffness J l s (T) are different from zero, as is the case at 
low temperatures, the average stiffness will be centered 
around the center of the Gaussian distribution (11), so 
that it will coincide with J (T). However, by approach- 
ing Tbkt defined by the average Jo(T) not all the patches 
make the transition at the same temperature, so that the 
BKT jump is rounded and J av remains finite above the 
average Tbkt, in agreement with the experiments, see 
Fig. 4. We note that Eq. (10) implies an average of the 
imaginary part of the complex optical conductivity, since 
J s oc A -2 ex f72- The same mechanism applied to its real 
part leads to a broadening of the <7i peak at the tran- 
sition, as discussed in Refs. [16,18]. Thus, one expects 
that a smearing of the abrupt superfluid-density jump 
due to increased inhomogeneity is also accompanied by 
a broadening of the o\ peak, as we will discuss indeed in 
the next section in connection to the experimental data. 



IV. DATA AND DISCUSSION 

The sheet superfluid densities d/X 2 (T) of many films 
with different thicknesses and disorder are shown in Fig. 
5. Fig. 5(a) shows moderately disordered films with 
T c > 8K. Fig. 5(b) are films with T c less than 8K. 
The modified BKT theory developed in Sec. Ill, with 
vortex core energy a free parameter and consideration of 
inhomogeneities, fit the experimental data pretty well. 
Table I shows major fitting parameters for every sample 
studied here. Two energy scales, the vortex core energy 
[i and the SC gap A(0), normalized to the superfluid 
stiffness, are found to be correlated at high disorder. This 
is shown in Fig. 6. The sheet superfluid density d\~ 2 (T) 
is converted to the superfluid stiffness J s by means of 
Eq. (2). By using /i 2 /4e 2 ^ofcs = 6.2 x 10~ 3 Km we can 
express J s in K as: 



J S [K] =0.62 



d[A] 
X 2 [fim 2 



(12) 



Qualitatively, all films show a deviation from BCS the- 
ory fit earlier than what 2D XY model predicts, which 
is the intersection with the dashed line in Fig. 5. This 
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FIG. 5: Experimental data (solid black curves) of various sam- 
ples [(a) for moderately disordered films and (b) for highly dis- 
ordered films] are fitted by BCS dirty limit theory well (blue 
dash-dotted curves) until deviations occur. Black dashed line 
gives the prediction on where BKT transition should occur 
given by 2D XY model. The deviations are fitted by the pro- 
cedure mentioned in Sec. Ill (red dashed curves). The fitting 
parameters values are reported in Table I. 



means that the early appearance of the BKT downturn, 
first shown in moderately disordered Nb 32 and NbN 5 
films, is a common characteristic of conventional 2D SC 
films. As we discussed in Sec. Ill, a small value of 
the vortex core energy is responsible for the fact that 
1/\ 2 {T BKT ) < 1/\ 2 bcs {T B kt), so that also the T B kt 
in the perfectly homogeneous case would occur before 
than the 2D-XY model prediction. However, the pres- 
ence of inhomogeneity smears out the sharp BKT drop 
and gives a finite width to the transition, as evidenced 
by small peaks in o\ near T c . These peaks get wider 
as disorder increases, consistent with the increase in in- 
homogeneity observed in tunneling 20-22 , and also with 
the increase of the superfluid-stiffness distribution width 
8/ J obtained by the BKT fit (see Table I). 

Quantitative results are shown in Table I. We are 



TABLE I: Experimental values of nominal thickness d and 
sheet superfluid density d/A 2 (0), which are transferred to en- 
ergy scale Js(0) via Eq. (2), along with the best fit parame- 
ters. Here the BCS transition temperature Tbcs and the su- 
perconducting gap A(0) are obtained from the BCS fit. The 
vortex-core energy /j, and the degree of inhomogeneity 8, both 
normalized by J 3 (0), are from best BKT fit. The tempera- 
ture Tbkt corresponds to the transition temperature of the 
median Jo(T) of the Gaussian distribution, see also Fig. (4). 



Film ID 


d 


d/A 2 (0) 


.7.(0) 


Tbcs 


Tbkt 


H/Js 


A/J 3 


S/Js\ 




(nm) 


(mm^ 1 ) 


(K) 


(K) 


(K) 








S56 


5.5 


59.4 


368.28 


13.7 


13.23 


0.6 


0.08 


0.007 


S87 


3.2 


24.96 


154.75 


12.23 


11.28 


0.54 


0.17 


0.013 


S84 


2.72 


15.2 


94.27 


11.1 


10.1 


0.7 


0.253 


0.02 


S88 


2.24 


8.31 


51.52 


9.7 


8.59 


0.83 


0.461 


0.02 


S194 


2.28 


5.68 


35.2 


9.3 


8.23 


0.75 


0.739 


0.04 


S196v2 


2.16 


4.02 


24.9 


7.7 


7.19 


1.65 


0.897 


0.05 


S197 


2.1 


3.93 


24.34 


7.95 


6.9 


1.25 


0.947 


0.045 


S196vl 


2.16 


3.82 


23.69 


8.68 


7.45 


1.55 


1.063 


0.045 


S193vl 


2.04 


3.41 


21.12 


7.48 


6.94 


1.7 


1.027 


0.042 


S193v4 


2.04 


2.69 


16.7 


7.1 


6.53 


1.9 


1.233 


0.055 


S193v5 


2.04 


2.58 


16 


6.78 


6.24 


1.93 


1.229 


0.062 


S193v2 


2.04 


2.04 


12.65 


5.87 


5.42 


2.2 


1.346 


0.07 


S193v3 


2.04 


1.41 


8.72 


5.15 


4.64 


2.25 


1.713 


0.085 


S89 


1.58 


1.20 


7.41 


4.2 


3.38 


1.8 


1.474 


0.08 



able to tune the disorder so that the superfluid stiffness 
changes by a factor of 50 and T c changes by a factor 
of 3. Our most disordered film has a T c ~ AK and 
J s (0) ~ 7AK compared to T c ~ 8K and J s (0) - 60K 
for the most disordered film in a previous study. 5 Tbkt 
is very close to T c at low disorder (3% difference) and 
becomes more separated as disorder is increased. For the 
most disordered film S89, this difference is as large as 
20%, showing the separation of two energy scales Tbkt 
and T c . Interesting, the energy scale for sheet superfluid 
density J s (0) is getting close to the scale of T c and may 
become the limiting factor for T c . Three fitting parame- 
ters, /i, A(0), 8, normalized by J s (0) are listed in Table 
I. 

One of our major findings is shown in Fig. 6. With 
J s (0) characterizing the disorder level of the films, we 
found that the two energy ratios A(0)/J s (0) and n/J s (0) 
are highly correlated at high disorder. This means that 
near the quantum critical point, the vortex-core energy 
/i, an important energy scale in BKT transition, does 
not scale with the superfluid stiffness, as given by 2D 
XY model. Instead, it scales with the superconducting 
gap, which is the pairing strength of the Cooper pairs. 

There are two different contexts our results can be put 
into. First, the observation of robust BKT transition 
is consistent with previous results on moderately disor- 
dered films 5 and extends it to highly disordered films on 
the verge of the 2D superconductor-insulator transition 
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FIG. 6: Evolution of the vortex-core energy (blue circles) and 
SC gap (red squares), normalized by the superfluid stiffness 
J s (0), with J s (0) in our NbN films. Open symbols are data 
from Mondal et al. Dashed lines are guides to the eye. 



(SIT). The vortex core energy is also shown to scale with 
the energy gap near the QCP. Our data on more than ten 
films firmly confirmed the observation of a previous study 
on three films, that can be understood theoretically 5 in 
terms of the increasing separation between the energy 
scales associated to pairing and phase coherence induced 
by disorder. 

Second, the robustness of BKT transition observed 
here is in direct contrast with similar superfluid den- 
sity studies on deeply underdoped layered cuprates 13-15 
near the QCP, where no downturns are observed at 
all. In these deeply underdoped quasi-2D cuprates, 
superfluid density goes linearly with the temperature 
in almost all the compounds studied, including both 
YBa 2 Cu307_.,; and Bi 2 Sr2CaCu208 +a; , both crystals and 
films. 

What is the big difference between the two systems? 
First of all, we should recall that the possibility to ob- 
serve a BKT transition in thick cuprate films relies on the 
interplay between different length and energy scales. In 
cuprates, there are three length scales in the c-axis: the 
thickness dfu m , the neighboring superconducting CuOi 
bilayer distance dcu0 2 an d the c-axis coherence length £ c . 
dfum is typically several hundred nanometers for "thick" 
films and much larger for crystals. The distance between 
neighboring C11O2 bilayers is about 12 A in YBCO and 15 
A in Bi-2212. Only the coherence length £ c has a temper- 
ature or doping dependence. Near well-studied optimal 
doping, £ c ~ 2A , which is much less than dcu0 2 - We 
have: 



£,c < dcuo 2 < d. 



film 



(quasi — 2D) 



(13) 



That is why cuprates are generally considered to be 
quasi-2D, and one would generically expect a BKT tran- 
sition for each isolated layer. This means that the tern- 



perature where the universal jump Eq. (9) should occur 
is compared to the supcrfluid stiffness of a single bilayer, 
i.e. d is replaced by dc u o 2 m Eq. (2). Nonetheless, lay- 
ers are not completely independent, since the phase in 
neighboring layers is coupled by a (weak) Josephson-likc 
coupling J±. Once more, when J± is much smaller than 
in-plane stiffness one would expect BKT-like behavior. 
However, this is only true when the vortex-core energy 
is of the value expected in the XY model: indeed, it has 
been shown 17 that for larger values of the vortex-core 
energy /i, the BKT transition can occur at temperatures 
larger than expected within the XY model, since inter- 
layer coupling predominates over vortex unbinding on a 
wider range of temperatures. The observation of sharp 
BKT downturns in optimally-doped Bi 2 Sr 2 CaCu208 +2 ; 15 
suggests that in this material not only is anisotropy very 
large (i.e. J±_ is very small), but also the vortex-core 
energy must be of the order of the XY-model value, al- 
lowing for a BKT transition controlled by the stiffness of 
each isolated bilayer. 

When doping is lowered both £ c and /i might grow. 
From one side, if £ c remains small, an increase of the 
ratio li/ J s analogous to the one reported above for con- 
ventional films could by itself move the Tbkt to higher 
temperatures, making also the jump barely visible. Such 
an increase of li/J s has been indeed inferred by a the- 
oretical analysis similar to the one discussed above in 
ultra-thin YBa2Cu307_ :c films. 24 In the case of conven- 
tional superconducting films the increase of /z/J s (0) can 
be understood as an effect of the increasing separation 
between the pairing energy scale and the phase coher- 
ence due to disorder, which can be also be responsible for 
the pseudogap observed by STM in this material. 19,20 A 
similar analysis in the context of cuprates would be very 
interesting, since it could shed new light on the effect of 
disorder on the underdoped regime of these materials as 
well. 

On the other hand, if £ c exceeds the dc u o 2 an d stays 
less than dfu m , 



1 CuQ 2 



<L< dju m (3D), 



(14) 



then cuprates become more three dimensional. A BKT 
paradigm then does not apply and the characteristic 
BKT jump disappears 33 . Eventually, when £ c exceeds 
even the film thickness 2D behavior could be recovered 
again, but now Tbkt should correspond to the sheet su- 
pcrfluid density of the whole film, which is very large for 
thick films. Thus, the BKT jump would become prac- 
tically indistinguishable from the transition temperature 
due to other thermal excitations (as quasiparticlc or lon- 
gitudinal phase fluctuations). The in-plane coherence 
length £ is related to upper critical field H C 2 by: 



H c2 = $ /2tt£ 2 



(15) 



where <J>o is flux quantum. Therefore a large coherence 
length corresponds to a relatively small H C 2 for deeply 
underdoped cuprates. In this view, the lack of BKT sig- 
natures in supcrfluid-density data can then support the 



idea that H c i drops and goes to zero near the underdoped 
side of superconductor-insulator transition. Of course, 
both the coherence length and the upper critical field 
are different along in-plane and out-of-plane directions, 
but we assume that this anisotropy is temperature and 
doping independent. This might provide some evidence 
on the recent debate 34,35 about how H c2 behaves on un- 
derdoped cuprates. At the QCP the coherence length di- 
verges for deeply underdoped cuprates and prevents BKT 
transition from occurring. 

The case of ultrathin conventional NbN superconduct- 
ing films is simpler. £ is isotropic and several tens of 
angstroms, to 2D behavior is always controlled by the 
thickness for ultrathin films: 



i>& 



film 



(2D), 



(16) 



Therefore they are always in the 2D limit and BKT tran- 
sition is always expected. Moreover, H c i measurements 
support the notion that £ increases when the QCP is 
approached, 36 so that even at strong disorder films re- 
main always in the 2D limit. Another way to think about 
the difference is, if we can reduce the thickness of cuprate 
films to a few unit cell so it is 2D by construction, then 
these films are similar to ultrathin conventional super- 
conducting films and we should be able to recover the 
BKT downturn. That is what we indeed see in ultrathin 
Y 1 _ a; Ca 2: Ba2Cu307_5 films. 27 In this case, a BKT-like 
downturn is observed and it is robust down to the lowest 
level of doping. 



V. CONCLUSIONS 

Robust Berezinskii-Kostcrlitz-Thouless transitions, ev- 
ident by sharp downturns of supcrfluid densities near 
T c , are observed for ultrathin NbN films close to 
superconductor-insulator transition. They occur earlier 
than what 2-D XY model predicts and this is attributed 
to a relative small vortex-core energy. We observe that 
the vortex-core energy scales with the superconducting 
gap instead of the supcrfluid stiffness near the quantum 
critical point. Once included this effect, the BKT transi- 
tion survives up to strong disorder, even though the sharp 
supcrfluid-density downturn observed in cleaner samples 
gets partly smeared out by the disorder-induced inhomo- 
geneity of the system. The robustness of BKT transition 
is in direct contrast to similar studies on severely un- 
derdoped layered cuprates, which show no critical ther- 
mal fluctuations. This difference is attributed to the ef- 
fect of a larger vortex-core energy or coherence length 
in deeply underdoped cuprates. Further investigation of 
both these mechanisms could shed new light on the na- 
ture of the superconductor-insulator transition in these 
unconventional superconductors. 
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